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Abstract 

We extend the canonical formalism for the motion of A^-particles in lineal gravity 
to include charges. Under suitable coordinate conditions and boundary conditions the 
Hamiltonian is defined as the spatial integral of the second derivative of the dilaton field 
which is given as a solution to the constraint equations. For a system of two particles 
the determining equation of the Hamiltonian (a kind of transcendental equation) is 
derived from the matching conditions for the dilaton field at the particles' position. 
The canonical equations of motion are derived from this determining equation. 
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For the equal mass case the canonical equations in terms of the proper time can 
be exactly solved in terms of hyperbolic and/or trigonometric functions. In electrody- 
namics with zero cosmological constant the trajectories for repulsive charges exhibit 
not only bounded motion but also a countably infinite series of unbounded motions 
for a fixed value of the total energy Hq, while for attractive charges the trajectories 
are simply periodic. When the cosmological constant A is introduced, the motion for 
a given A and Hq is classified in terms of the charge-momentum diagram from which 
we can predict what type of the motion is realized for a given charge. In general the 
cosmological constant acts on the particles as a repulsive (A > 0) or an attractive 
(A < 0) potential. But for a certain range of A < 0, small q and small mass the 
trajectory shows a double peak structure due to an interplay among the induced mo- 
mentum dependent A potential, the gravitational attraction and the relativistic effect. 
For A > 0, depending on the value of charge, only bounded motion or the infinite series 
of unbounded motion or both are realized. 

Since in this theory the charge of each particle appears in the form eie2 in the 
determining equation, the static balance condition in 1-1-1 dimensions turns out to be 
identical with the condition in Newtonian theory. We generalize this condition to non- 
zero momenta, obtaining the first exact solution to the static balance problem that 
does not obey the Majumdar-Papapetrou condition. 
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1 INTRODUCTION 



The problem of determining the motion of a system of particles mutually interacting 
through specified forces is one that has attracted attention since the dawn of physics. We 
continue in this paper our previous explorations of this problem in two spacetime dimensions 
when the specified interactions are both gravitational and electromagnetic. 

Our work here represents the first exact solution to include both interactions in a rela- 
tivistic framework. Although an exact solution is known in three spatial dimensions for pure 
Newtonian gravity in the N = 2 case, dissipation of energy in the form of gravitational ra- 
diation has obstructed progress toward obtaining exact solutions for the motion of bodies 
in the general theory of relativity, even when N = 2. However by reducing the number 
of spatial dimensions this obstruction disappears, at least in the vacuum. Apart from the 
absence of gravitational radiation, most (if not all) of the remaining conceptual features of 
relativistic gravity are retained, and so lower dimensional theories of gravity offer the hope 
of garnering insight into the nature of classical and quantum gravitation in a wide variety 
of physical situations. 

For these reasons we have been investigating the A^-body problem in two dimensional 
gravity for the past 3 years. We have chosen to work with a 2D theory that models 4D 
general relativity in that it sets the Ricci scalar equal to the trace of the stress-energy 
of prescribed matter fields and sources. Hence matter governs the evolution of spacetime 
curvature which reciprocally governs the evolution of matter This theory (sometimes 
referred to as i? = T theory) has a consistent Newtonian limit a problematic limit in a 
generic (1 + l)-dimensional theory of gravity theory [Q. When the stress-energy is that of a 
cosmological constant, the theory reduces to Jackiw-Teitelboim (JT) theory p. 

Working in the canonical formalism Q , we have so far been able to obtain exact solutions 
to the two body problem in the absence and presence P of a cosmological constant. 
In this paper we extend these considerations to include charged bodies. Specifically, we 
formulate the charged A^-body problem in relativistic gravity by taking the matter action 
to be that of A^ charged point-particles minimally coupled to gravity and electromagnetism. 
We extend the previous canonical formalism we developed for this action in i? = T lineal 
gravity to include this case. When N = 2 we obtain exact solutions for the motion of 
two bodies of unequal (and equal) charge and mass. In the slow motion, weak field limit 
the Hamiltonian we obtain coincides with that of Newtonian gravity with electromagnetism 
in (1 + 1) dimensions. We are also able to extend our solutions to include a cosmological 
constant A, so that in the limit where all bodies are massless and neutral, spacetime has 
constant curvature (ie the JT theory is obtained). Our solution is the first non-static exact 
solution to the charged 2-body problem in any relativistic theory of gravity. 

Our exact solution to the N = 2 case can be formulated in several ways. We derive an 
exact solution for the Hamiltonian as a function of the proper separation and the centre-of- 
inertia momentum of the charged bodies. We are also able to construct a solution in which 
the proper separation of the two charged point masses is given as a function of their mutual 
proper time in the equal mass case. If the masses are not equal our exact solution is given 
in terms of a time-coordinate that is not the proper time. 

A scalar (dilaton) field must be included in the action since the Einstein action is a 
topological invariant in (1+1) dimensions. Canonically reducing the action, we find that the 
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Hamiltonian is given in terms of a spatial integral of the second derivative of the dilaton field, 
which is a function of the canonical variables of the particles (coordinates and momenta) 
and is determined from the constraint equations. The matching conditions of the solution to 
the constraint equations yield an equation which determines the Hamiltonian in terms of the 
remaining degrees of freedom of the system when N = 2 . We refer to this transcendental 
equation as the determining equation, since it allows one to determine the Hamiltonian in 
terms of the centre of inertia momentum and proper separation of the bodies. From this 
Hamiltonian we can derive the canonical equations of motion. In the equal mass case we find 
that the separation and momentum are given by differential equations in terms of the proper 
time, and can be exactly solved in terms of hyperbolic and/or trigonometric functions. 

Several different types of motion are expected in the 2 body case, depending upon the 
signs and magnitudes of the masses, charges, energy and other parameters (e.g. gravitational 
coupling constant, cosmological constant). If the charges are of opposite sign the particles 
will remain bounded, whereas if they are of the same sign either bounded or unbounded 
motion can occur for the same value of the total energy. For a given set of parameters there 
is in this case a countably infinite series of unbounded motions labelled by an integer n. 
A balance condition exists between the bounded and the unbounded cases, and reduces to 
the expected (Newtonian-like) static balance condition in the absence of particle momenta. 
We shall analyze these various states of motion, and discuss the transitions which occur 
between them. A cosmological constant can qualitatively change the motion, rendering 
bound states unbound and vice-versa. We find several surprising situations, including the 
diverging separation of the two bodies at finite proper time in the repulsive case, even for 
A = 0. In the A < case the motion shows a double maximum behavior for a certain range 
of the parameters, which is a characteristic effect of the induced momentum dependent A 
potential. For classification of the motion we utilize a charge-momentum diagram from which 
we can easily predict what type of the motion is realized for a given charge. 

In the unequal mass case the proper time is no longer the same for the two particles, and 
a more careful analysis is necessary in order to describe the motion. We obtain phase space 
trajectories from the determining equation and explicit solutions for the proper separation 
in terms of a transformed time coordinate which reduces to the mutual proper time in the 
case of equal mass. 

In Sec. II we describe the outline of the canonical reduction of the theory when the 
charges are included and define the Hamiltonian for the N- body system. In Sec.III we 
solve the constraint equations for the two-body case and get the determining equation of 
the Hamiltonian, from which the canonical equations of motion are explicitly derived. For 
the case of equal masses and arbitrary charges the explicit exact solutions to the canonical 
equations are given in Sec. IV. By using these exact solutions we analyze in Sec.V the motion 
for A = separately for four cases: attractive charges, small repulsive charges, large repulsive 
charges and H < 2m. We analyze the motion of equal masses for A in Sec. VI for four 
possible combinations of the signs of A and the charges, where we also develop the plots 
of phase space trajectories and the analysis of the explicit solutions in terms of the proper 
time, based on a classification given by a charge-momentum diagram. We treat the unequal 
mass case in Sec.VII. In Sec.VIII we investigate the static balance problem by using both 
the canonical equation and the determining equation. Sec.IX contains concluding remarks 
and directions for further work. The solution of the metric tensor, a linear approximation of 
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the exact solutions and the causal relationships between particles in unbounded motion are 
given in Appendices. 



2 CANONICALLY REDUCED HAMILTONIAN OF 
iV-CHARGED PARTICLES 

Our derviation of the canonically reduced Hamiltonian for charged particles is parallel to 
that given in the uncharged case 0, ^. Here we briefly review this work, highlighting those 
aspects that are peculiar to the charged case. 

The action integral for gravitational and electromagnetic fields coupled with charged 
point masses is 
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where \l/ is the dilaton field, and J-'^'^ are the vector potential and the field strength, g^^^, 
and g are the metric and its determinant, R is the Ricci scalar, and and are the charge 
and the proper time of a-th particle, respectively, with k = SttG/c^. The symbol denotes 
the covariant derivative associated with gfj_^. 



The field equations derived from the action (|T|) are 

R - g^'V.V,^ = , 

= E ea / dTa^6\x - zM) , 
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dTa 



dz^ 

dTr, 



/J. ) 



dz'^ 

dTn 



(2) 
(3) 

(4) 
(5) 
(6) 



where the stress-energy due to the point masses and the electric field is 



f 1 dz^ dz^ If 1 1 

T^u = E"^W dTa^=g^.g,p-^-^6\x - zM) + -. — r T^^T,fig''<' - -g^.uT^fiT''^ , 

a J \'-9 dTa dTa {-g) I 4 J 

(7) 

recalling that in (1+1) dimensions no magnetic component of the field exists. Eq.(^ guar- 
antees the conservation of T^^- Inserting the trace of Eq.(|^) into Eq.(0) yields 



R-A = kT^",, 



Eqs. (H) , @ , (I) and (H) form a closed sytem of equations for gravity, electromagnetism, and 
matter. 

On transforming the action (|l|) to the canonical form we first rewrite the electromagnetic 
and the particle sectors in a form appropriate to the first-order formalism P[ . 



lE+P 
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where vr^ is essentially the four velocity of the particle and is a Lagrange multiplier. The 
variations with respect to Ay, and T^^ lead to Eqs.(Q) and (^), respectively, and those with 
respect to tTq^, and A'^ lead to 



dTa 

dn, 



K< = , 
1 



dz'i 



+ 7^K'^a\T^av9^'^y{Za) = {A^,f,{Za) - Af,^,y{Za)} 



T^ai,'^av9 + "^a = • 

This set of three equations is equivalent to Eq.(^ when A'^ = l/rria is chosen. 
By performing the integration over the parameter and setting 



(10) 

(11) 
(12) 



(^, Za) 



-^,A) 



Zl^ = (t,Zn] , IT, 

A, = (- 

\ _ \' I 



E = r 



01 



(13) 
(14) 



the action (M) becomes 



lE+P 



X - Zait)) 



jd^x\A---V^E^ + ^l^--^eJi 
+ E + eaA)5(a; - Zait)) + E {^Tao - ^A(7r,^7r,,(7'^" + m^)| 5(x - 2;,(t)) 



Varying the Lagrange multipliers if and A^ yields the constraints 

dE 



dx 



^ea5{x - Za{t)) 



The solution to (ITB) is 



^^a^l^^au9^'' + m„ = . 



(16) 
(17) 



In (1+1) dimensions the electric field has no independent degrees of freedom. When the 
charged particles move within a finite region [\x\ < L), the electric field in an outer region 



5 



> L) is E = i^Z^a^a- For a system of zero total charge the electric field vanishes in 
the outer region. As the solution to (^) we choose 

VTaO = ^ { + V(^?°'7rJ2_ ^00(^11^2 + ^2)| . (ig) 

After the constraints are eliminated, the action (|T^) is 



l-E+P 



dz 

j:rr^^6{x - z^{t)) (20) 

a 9 ^ J ^ 

From this expression we know that TTa is the conjugate momentum to Za and hereafter we 
use the notation 

Note that the variations of the action (|T5]) with respect to A and E lead to 

^ + i:e.^'5(x-z„(i))=0, (21) 

Eg . (1211) is automatically satisfied for the solution (|18D to the constraint equation, in contrast 
to the (3+1) dimensional case where the corresponding equation is a true dynamical equation. 
Upon inserting the solution ([T8|) into the action (|T5|) all terms related to the components of 

cancel. Hence we need no longer consider A^. Its solution, if desired, is straightforwardly 
obtained by solving ([16D after fixing the gauge and obtaining the full solution of the metric. 

Consider next the transformation of the gravity sector to canonical form. We decompose 
the scalar curvature in terms of the extrinsic curvature K via 

^gR = -2d^{^K) + 2d^[{N^K - d^N^)/^] , (23) 

where the metric is ^ 

ds"^ = -N^dt^ + -i[dx + ^dt^ , (24) 

with K = {2No^)-\2d,m-^-'N,da-do7), so that ^ = gn,No = (-^°°)-^/' and N, = gm- 
We then rewrite the gravity sector in first-order form . We find that the action ([^) becomes 



(25) 



J = jd'^X \^PaZa5{x - Zait)) + 717 + + N^R'' + N^R"^^ , 

where tt and 11 are conjugate momenta to 7 and \E', respectively, and 

a \ I 

R^ = -vr- in*' + 27r' + y ^(5(x-;2,(t)) , (26) 

7 7 a 1 
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with the symbols ( ' ) and ( ' ) denoting do and di, respectively. 
From the action we obtain the set of equations 



TT 



Nol^iVi-' - ^-n + + T^iE' - -) 

-E ^—6{x-z^{t)) 



7 7 a 7 J 2fi;^7 7 



7-Aro(2/€V777r-2/€7^) + iVi^-2Ar( = , (28) 

7 

= , (29) 

i?^ = , (30) 

n + 9i(--iVin + -^iVo^' + -^K) = o , (3i) 

7 2k^ 

^ + Aro(2fi:v^vr) - Ari(-^') = , (32) 

7 



Pa + —+^a ; , — o 1^ 

•^^a V 7 2 / ^ + m2 ^ ^^"^ "^^^ ^ 



+ '^'^ai + /''^W5;; = 0, (33) 

iVo , ' + — = . (34) 

7 

In the equations ( ^3l) and (|3^), all metric components (A'o, A''i, 7) are evaluated at the point 

X = Za and 

df _ df{x) 



dza dx 

The quantities Nq and Ni are Lagrange multipliers which yield the constraint equations ( p9|) 
and (^). The above set of equations can be proved to be equivalent to the set of equations 
(I), (I) and (i. 

To proceed to the canonical reduction of the action (P3| ) we have to eliminate the redun- 
dant variables by utilizing the constrant equations to fix the coordinate conditions. Noticing 
that the only linear terms in the constraint equations ( |29| ) and ( |30D are {^!' / ^/^^ and vr', 
respectively, and the equations may be solved for these quantities, we can transform the total 
generator obtained from the end point variation into an approriate form to fix the coordinate 
conditions. Generalizing the procedure described in our previous papers 0, |[ for the case 
of A = Cq = 0, we find that we can consistently choose the coordinate conditions 

7 = 1 and H = . (35) 
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Eliminating the constraints, the action (|25|) reduces to 



(36) 



where the reduced Hamiltonian for the system of particles is defined by 



H = I dxH = -- / . (37) 



K 



Here \l/ is a function of Za and pa and is determined by solving the constraints which under 
the coordinate conditions (B^) become 



Avl/ - -im'f + k^t:^ + -{kE^ - A) + + ^2^(3, _ ^ J = Q , (38) 

27r' + ^p,5(x-;2j = 0. (39) 

a 

The consistency of this canonical reduction is proved in an analogous way to the case of 
A = Ca = : namely the canonical equations of motion derived from the reduced Hamiltonian 
(|37D are identical with the equations (^) and (Q) [|, |]. 

3 SOLUTION TO THE CONSTRAINT EQUATIONS 
AND THE HAMILTONIAN FOR A SYSTEM OF 
TWO PARTICLES 

The standard approach for investigating the dynamics of particles is to get first an explicit 
expression of the Hamiltonian and to derive the equations of motion, from which the solution 
of trajectories are obtained. In this section we show how to derive the Hamiltonian from 
the solution to the constraint equations (53) and (BO) and get the explicit Hamiltonian 



for a system of two charged particles. Since the electric field appears in the combination 
[E"^ — A/k) in all equations we set 

V{x) = E^~C and A^ = A - /tC (40) 

with C = ^ij^a^a)^- Thus V{x) vanishes in the outer region and Ag is an effective cosmo- 
logical constant, which includes the contribution from the electric field. This latter situation 
arises from the well-known fact that in (1+1) dimensions the electromagnetic field strength 
is a 2-form, and so in compact spatial regions it contributes to the stress-energy tensor in the 
same manner as a cosmological constant, analogous to the way in which a 4- form behaves 



in 3-1-1 dimensions |]T0|- We shall later see that when Ag vanishes we get the Hamiltonian 
which leads, in the limit k — 0, to the correct special-relativistic electrodynamics in (1+1) 
flat space-time. 



We express equations (|38|) and (|3^) as 



Avl/ = i (M/')' - ixf - li^^V - A,) - KY.^pl + ml6{x - z^) , (41) 
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^X = -^EPa5(a;-Za), (42) 



where we set x' = ^- Rewriting (^Tj) as 



(1 + -) Av[/ = - A(vI/2 - Ak\^) -^{kV- K) + k\/\x -'^T.^pI + - ^a) , (43) 

yields, upon insertion into the Right-hand side (RHS) of ( P7| ) 



4^ 1 + ivl;(,J + 2 4^ 1 + ivl/(z,) + 2 i ^ + ivl/(xj 

-— / dx ir^A f^^ - 4/«\2 + 2Aea;2) . (44) 

8k J l + i^(x) V A e y V ; 

an expression which can also be obtained repeated iteration of the insertion of the RHS of 
(|Tp into the RHS of (^. 
Defining (p by 

^ = -4 In 101, (45) 



the constraints ( pj) and (|4^) for a two-particle system become 



A0-i{«:2(x')' + V-iAe}0 = ^{Vp? + mf0(^i)5(x-^i) 

+ ^Jpl + ml (I)iz2)6{x - Z2)] , (46) 



Ax = - 2:1) + P2S{x -Z2)} . (47) 



The general solution to (^Tf ) is 

X = -J {pi I a; - zi I +P2 I X - Z2 1} - e^a; + eC^ . (48) 

The factor e (e^ = 1) has been introduced in the constants X and so that the T-inversion 
(time reversal) properties of x are explicitly manifest. By definition, e changes sign under 
time reversal and so, therefore, does x- 

Consider first the case Z2 < Zi, for which we may divide space into three regions: Zi < x 
((+) region), Z2 < x < Zi ((0) region) and x < Z2 ((-) region). In each region, V and x' are 
constant: 

{0 (+) region, 

-ieies (0) region, (49) 

(-) region, 

-eX ~ jipi+ P2) (+) region, 

X^ = { -tX + i(pi -P2) (0) region, (50) 

-eX +4(^1+^2) (-) region . 
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General solutions to the homogeneous equation A0 — j (x')^ + — ^Aej = in each 
region are 



where 



k2 (x + |(pi+p2; 



/t2 (X- f(pi -p2)j -feiea 

(X-|(pi+p2))'-iAe 



(+) region , 
(0) region , 
(-) region . 



(52) 



For these solutions to be the actual solutions to Eq.(^6|) with delta function source terms, 
they must satisfy the following matching conditions a,t x = zi, Z2: 



M^i) = 0(^1) > 
00(^2) = 0(^2) , 

00(^2) - 0'_(2;2) = 7VP2 + ™ 



^1 



10(2:2) 



(53a) 
(53b) 

(53c) 
(53d) 



Since the magnitudes of both \l/ and x increase with increasing it is necessary to 
impose a boundary condition which guarantees that the surface terms which arise in trans- 
forming the action vanish and simultaneously preserves the finiteness of the Hamiltonian. 
From the iterative expression (|^) we know that we may choose the boundary condition 

\1>^ - 4k^x^ + 2AeX^ = C±x for (+) and (-) regions (54) 

with C± being constants to be determined. 

The above matching conditions accompanied by the boundary condition (|5^) determine 
the solution ( and also all coefficients) completely. The process of the calculation is quite 
analogous to the procedure in the previous paper , and we shall omit the details here. The 
compact expression of the solution is 

( Q-\{KoiZ^-Ko2Z2)+\K+(x-Zi) 



\Mi. 

— 6-3(^01^1-^02^2) |(iriA^i)i/2g-ii^o(x-.i) ^ (i^2A<2)^/^e^^o(^-^2)| , (55) 



4Ko 



e 



- i (KoiZi-Ko2Z2)-^K- {X-Z2) 



where 



Ki = 2Ko + 2K^ - nJpi + mi 



K2 = 2Ko + 2K+- nJpl + ml 
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i^oi = Ko-K+ + ^pi, (56) 



and among these quantites there exists one relation 

KiK2 = M1M2 e^o("^-"^) , (57) 

which we refer to as the determining equation of X. 
The Hamiltonian (pTl) becomes 

If 4 
H = / dx/\^ = - 

L J —00 

(5^ 



K J K 

2{K+ + K„) 



AC 



Once the solution X to ([57| ) is obtained, the Hamiltonian is explicitly determined. Conse- 
quently ( |57D is the determining equation of the Hamiltonian. 

Repeating the analysis for zi < Z2 yields a similar solution with pi — >■ —pi and so the 
full solution is obtained by replacing pi and zi — Z2 by pi = Pi sgn(2;i — Z2) and \zi — Z2\, 
respectively. The determining equation (|57|) of the Hamiltonian becomes 

K1K2 = M1M2 e^"l"^-"^l . (59) 

or more explicitly 

Uk^ + [k^JpI + ml - 2K_^][K^pl + ml - 2KJ\^ tanh (^^oki - ^^al) 



= -2K^ ([k^pI + mi- 2K^] + [n^pl + mi - 2KJ\^ (60) 

where the momentum pi is replaced by pi. 

For the expression ( |5BD to have a definite meaning as the Hamiltonian, K± must both 
be real, with positive sum. This imposes a restriction on X corresponding to a value of the 
cosmological constant Ag. However Kq is not necessarily real. If 6162 takes a large positive 
value (strong electromagnetic repulsion), Kq may be imaginary. In this case we need to 
reconsider the above analysis, because in the (0) region the soluton to the equation ( ^6|) 
becomes ^ ^ 

(j)Q{x) = As sin-^oa; + cos-i^oa: , (61) 

where 



'^e,e2 + lAe-^'{x-'-{p,-P2)) . (62) 
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Under the same matching conditions (|53a[ - |53d| ) and the boundary condition (|5^) we get, 
instead of (|5UD, a new determining equation of the Hamiltonian 

•1 



- 2KJ\ ) tan ( -Kq\zi 



Z2\ 



2Ko [kJpI + mf - 2K+] + [^Jpl + ml - 2K 



(63) 



Actually, this is just the equation derived from ( |5D[ ) by formally replacing Kq with iKo. The 
solution of (j) for imaginary Kq is also identical with that obtained from (|55D by the same 
replacement. We can therefore use equation ( ^9]) for all values of Ko; it is a transcendental 
equation which determines H in terms of the momenta and positions of the particles. We 
have previously shown that in the case of zero cosmological constant and no charges the 
solution for H is expressed in terms of the Lambert W function. More generally, with Ag 
and Cq all nonzero, a solution to (|59|) for H cannot be explicitly expressed in terms of known 
fuctions. However it can be obtained in successive approximations in the parameters c~^, k, 
etc. The examples will be shown later in the sections VI and VIII. 

Though in the general case the Hamiltonian cannot be expressed explicitly in terms 
of known functions, the canonical equations of motion can be exactly derived from the 
determining equation (|59|) by differentiating it with respect to the variables Zi and pi . For 
the variables pi and zi we have 

dH 2 fdK. 



Pi 



+ dK_ 
yi \ KoKiK, 





(64) 







+ 2^- 




\ dX 




-a 






dpi 



8 fY+ , KqKi 
Ml 



Pi 



Ipi + mf 



(65) 



where 



y+ 

Y_ 



K 



K 



K 



X + -{pi+P2) 

X - ^{pi -P2) 

X-^-{pi+P2) 



(66) 



and 



J 





(67) 
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Similarly for particle 2 the equations are 



P2 



Z2 



— e 













-) J ' 






Y+ Y\ K0K2 


=7r + 


■'A 


K_) M2 



P2 



/p2 + m 



(68) 
(69) 



These canonical equations guarantee the conservation of the Hamiltonian and the total 
momentum pi + P2- 



On the other hand the equations of motion (|33D and 
become under the coordinate conditions ( ^5]) 



derived from the action ( p5D 



Pa 



dZa 



pl + ml + 



dNi 
dza 



'P^ + 9 Z^^a^b^O^Fa - Zb\ 



Pa 



pl + m_ 



(70) 
(71) 



It is straightforward to verify that insertion of the solutions of the metric components 
given in Appendix A into ( [TOD and ([7T|) reproduces the canonical equations of motion (0), 
and (B^) where the partial derivatives at zl, z2 are defined by 



giVo,i 
dzi 



dNo,_ 



dx 



+ 



dx 



(72) 



Thus consistency between the geodesic equations and the canonical equations of motion is 
explicitly assured, while formal proof of consistency in the case of Ag = and Ca = can be 
easily generahzed 0]. 

The components of the metric are determined from the equations (]27|), (|28|), (|31| ) and 
(|32[) under the coordinate conditions (|35|). The derivation and the explicit solutions of the 
metric are given in Appendix A. With these solutions we can trace how the structure of 
space-time changes due to the motion of the two bodies. 



4 EXACT SOLUTIONS OF THE TRAJECTORIES 
FOR EQUAL MASSES AND ARBITRARY CHARGES 

In this section we consider a system of two particles with equal mass. Since the total 
momentum is conserved we can always choose the center of inertia (C.I.) system pi = —p2 = 
p. Corresponding to the sign of {^JJp'+'sAJk? — 2ep)^ — 86162/ k. — SAg/ the determining 
equations ( [60| ) and ( |63D become 

(J/ + B^) tanh (^^ Ja |r|) = 2JaB , (73) 

and 

iJ^^ - B^) tan {^^J^ |r|) = -2J^B , (74) 
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respectively, where 



2ep 



86162 8Ag 



K 



K 



2 ' 



^6162 8Ae 



IJJ2 



8Ae 



2ep 



B = H -2^p^ + m? . 
The equation (|73D is further divided into two types: 



tanh 



16 



MA 



(tanh-type A) 



or 



tanh ( y^^akI 



B 

^ , (tanh-type B) . 



(75) 



(76) 



(77) 



In the case of Ag = and 6a = the tanh-type B equation is excluded, because JkI B 
exceeds unity. When a cosmological constant and/or charge are introduced, this equation 
may have a solution in some restricted range of the parameters. 
Eq.([7^) is also divided into 



or 



tan ( -^Jk\r\ 



tan ( Y^^Arl 



B 
Ja 



, (tan-type A) 



^ , (tan-type B) . 



(7^ 



(79) 



For all four types of the determining equations the canonical equations of motion are 
identical: 



P 



52 



sgn(rj 



^0) 



2e Wl 



8Ae 



1 - 



c 



sgn(r) + 



P 



C ^/pP^+m? 



where 

C = 



8Ae 



'1 + 



8Ae 



Jl 



B^) r 



12) 



For given values of Ag and 6a, the equations (|73| ) or (|7^ describe the surface in (r, p, H) 
space of all allowed phase-space trajectories. Since if is a constant of motion, we can draw a 
phase space trajectory in (r, p) space by setting H = Hqih (^) or (|74|). This same trajectory 
can also be obtained directly from the solution r{t),p{t) to the canonical equations (^) and 
(|8lD by eliminating the time variable t. Numerical solutions to (^) and (^Tj) confirm this, 
but in r(t) and p{t) superficial singularities appear due to the zero points of C. 



14 



It is therefore preferable to describe the particles' trajectories in terms of some invariant 
parameter. The common proper time of each particle is the best candidate as seen in the 
starting action (|I]). From the metric components given in the Appendix A and the canonical 
equations (j7T|), the proper time is 



Tfl 

= deN,{z,f-—^^ (a = 1,2) . (83) 

For the equal mass case it is identical for particles 1 and 2 

dT = dn = dT2 = -^^=^dt , (84) 
from which the canonical equations (|80|) and (pl|) may be expressed in the form 



2 



dp nVp^T^jJl - B') 

-r = sg'^(^) ' (85) 

dr Ibm 

dr 2e fv^p^Tm^C , ^ ~J 

_|_ _ y sgn(rj . (8oj 



2 



dr m \ Jj 

Remarkably these equations can be solved exactly. First we solve Eq.(|85|) for p(r) and 
then obtain r(r), either by directly solving (|86D after substituting the solution for p or by 
solving ([761) -([79|) for r. For the r > region Eq. ([85|) leads to 



p dp kH 

kH 



dr 



This expression implies the condition 



l + |^>0 (88) 
which is satisfied for all Ag > 0. For negative Ag the motion is allowed provided H satisfies 



8Ae 

H>f-^. (89) 



We perfrom the integration of the Left-hand side (LHS) of (|87D in three cases, separately, 
depending on the value of A^. The solution p(r) is 

Pir) = f (fir) - ^) (90) 
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with 



fir) 



^ 1 



l+VTfi" 

771 I (7 



/7m{T-To) 



--J^ V— 7m tan 



-7m(T-To) 



-7m(T--To) 



7m > , 
7 m = , 

7m < , 



where 



7h 

7m 
^ = 



1 _L 8Ae 



7e = 1 + 



2eie2 



(1 + v/7^)(v/pg + m2 - epo) - ^ 7e 



with po being the initial momentum at r = tq. 
Similarly the solution in r < region is 



with 



1+v^7h" 



/7m(T-rg) 



7e+V-7m— Y 



-7m{T--To) 



-7m(T--To) 



where 



(X = (1 + ^)Nvl + m? + epo) 



m 



■7e 



7m > , 
7m = , 

7m < , 



(91) 



(92) 



(93) 



(94) 



(95) 



Corresponding to each type of the determining equations ([761)- ([79| ) the solution for r(r 
is obtained as follows, 



16 



tanh-type A: 



r(T) = 



16 tanh ^ 






7^1) 




\ 


/ (^^yK■2H2+8Ae-mK{f{T)- 


7m)^ 


2 

-8fteie2-8Ae 



-16 tanh ^ 



(^^K2j^2+8Ae-mK(/(r)--^)^ -8Keie2-8Ae 



tanh-type B: 



r r 



16 tanh ^ 


^Vi^2-ff2+8Ae-mK(/(T)- j^))^-8fteie2-8Ae 




«-)+7f?)l) 


\ 


/ (y K^H^+8Ae-mK{f{T)- 


\ 2 

-8/teie2-8Ae 



-16 tanh ^ 



y ^V«2if2-H8Ae-mK(/(T)--=i^) 


1 — 8Ke]^e2— 8Ae 






) 



2i?2+8Ae-mK(/(r) 



1 ^ 



8«;eie2— 8Ae 



tan-type A: 



r(T) 



16 tan ^ 



8Ae-|-8Keie2-(-\/K^w2+8X^ 



-l-n-TT 



-16 



tan ^ 



8Ae-f8Keie2-(\/K2i?2+8Ae-mK(/(r)--^) 

") 



S.\, ■ r<i 



v'^■-^r.f-^-:-8A- mK(/(T) - -i--) 



-|-n7r 



8Ae-f8Keie2- if2+8Ae -mK(/(r) - )) 



r > , 



r < , 



r > , 



r < , 



r > , 



r < , 



(96) 



(97) 



(98) 
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tan-type B: 



r(r) 



16 tan ^ 



8Ae+8fceie2-(VK2H2+8Ac-mK(/(T)-ji:j))^ 

«(ff--|/M+7f^|) 



y8Ae+8Keie2-(^A/K2/^2+8Ae-mK(/{T)-yi^)^ 



'16 



tan ^ 



+ 8Ae-mK(/{T)----!— ) 



T 

1 



8A. 



;+8Keie2- (^v^K2H2+8Ae-mK(/{T)-ji^) 



r > 



(99) 



r < 



5 ANALYSIS OF ELECTRODYNAMIC MOTION WITH 

Ae = 

Based on the exact solutions in the previous section we can analyze the dynamics of two- 
body problem. In this section we investigate electrodynamics in a space-time with Ag = 0. 
The solution p(r) in this case is given for r > by 



with 



and for r < by 



with 



dTTl I 1 



fe{r) 



p{t) 



Mr) 



H 

m7e 



pl + m'^ - epo 



1 



H 



Po + m'^ + epo - 



PQ + m^ + epo 



(100) 



(101) 



(102) 



(103) 



The relative distance r(r) is obtained from (76) - (|T9D by simply replacing /(r) and /(r) 
with /e(T) and /e(T), respectively. 
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Relative motion of two charged particles is classified by the signs and the magnitudes 
of the charges. Since the charges always appear as the product 6162 in all solutions, it is 
sufficient to analyze the attractive case by setting ei = —62 = q and the repulsive case by 
setting ei = 62 = g so that the charges have equal magnitude. Throughout this paper, in 
the numerical analysis, we set e = 1, k = 1 and rescale everything in units of m (effectively 
setting m = 1) for simplicity, except when otherwise stated. It should be remarked that for 
every phase space trajectory there exists a time-reversed trajectory with e = — 1, which is 
obtained by reversing the former with respect to the r-axis. 



5.1 The attractive case: ei = —62 = q 

When the electric force between charges is attractive, the value of (^H — mf{T) + j^'j 
8^6162 is always positive and two particles follow a bounded periodic motion described by 



tanh-type solution (|76|). The period is determined from the initial value of po = y {H/2y — w? 
at r = 0: 



nm^e \ (2 - ^e)H 

Although the above expression diverges when 7e = 2p^H ' ^^^^ situation is never realized in 
the attractive case, since 7e < 1 whereas 2p^H ^ ^■ 

In Figs.l and 2 we plot r(r) and the phase space trajectories , respectively, for two 
particles initially at r = for fixed = 1 and four different values of H. We see that as H 
increases, the phase space trajectory becomes more S-shaped, r(r) becomes steep and the 
amplitude increases, but the period approaches a constant value (16/Km7e) tanh~^[7e/(2 — 
7e)]. This S'-shaped deformation of the trajectory at higher energy is caused by the p- 
dependence of the gravitational potential and is common to the g = case described with 
the W function ||^. Figs. 3 and 4 show similar plots for fixed H = 3 and four different values 
of |g|. For large |g| (7e < 0) due to the attractive force between charges the phase space 
trajectory contracts toward the origin, and the period and the amplitude of r(r) become 
small. As |g| becomes small, 7e passes zero at |g| = and approaches 1 at |g| = 0. As 

we can see from the figures, no peculiar changes in either r(r) or the phase space trajectories 
occur except for a growing amplitude and period, a natural tendency due to the weakening 
attractive force. 

The distinct characteristics of these relativistic motions are more easily seen by comparing 
the exact trajectory with those of the motions in three approximations : 

(1) the non- relativistic motion described by the Hamiltonian 

H = 2m + ^ + ^\r\ + ^\r\, (105) 
m 2 4 

(2) the linear approximation in k, its Hamiltonian being 

+ \<f{^P^ + m2 - ep)r^ + , (106) 
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Figure 2: Phase space trajectories correponding to the r(T) plots in Fig.l. 
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Figure 3: The r(T) plots for Hq — 3 and four different values of \q 



22 



p 




Figure 4: Phase space trajectories correponding to the plots in Fig. 3. 
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(3) the K — > limit theory, namely, special-relativistic electrodynamics in (1 + 1) flat space- 
time described by the Hamiltonian 

2 

H = 2^p^ + m? + '^ \r\ . (107) 

In Fig. 5 the trajectories of the exact and the three approximate solutions under identical total 
energy Hq = ?> are drawn for q = 0.5, 1, 5 and 10. For small q = 0.5 both the exact solution 
(solid curve) and the hnear approximation (dashed curve) follow the S'-shaped trajectories, 
while the trajectories of the non-relativistic solution (dotted curve) and flat electrodynamics 
(dot-dashed curve) have symmetrical oval forms. As |g| increases all trajectories tend to 
coincide with the trajectory of flat electrodynamics (since the effect of gravity becomes 
relatively weak) though the initial value of po = \Jm{H — 2m) = 1 in the non-relativistic 
case is slightly different from the value of others (po = \/ {.H/2Y — m? = a/5/2). If the 

motion in the non- relativistic case starts from the same initial value po = the oval 

trajectory for large |g| region as shown in Fig.6 becomes larger than others, reflecting the 
difference between the relativistic and non- relativistic effects. 



5.2 The repulsive case: the small \q\ regime (ei = 62 = q, \q\ < qc) 

For the case where the charges are of the same sign, the electric force is repulsive and 
competes with the attractive gravitational force. Depending on the strengths of the two 
forces the solutions become either tanh-type and/or tan-type. The transition from tanh- 
type to tan-type is given by the zeroes of = ^ {H — 2epy — Sg^/ k, which lead to a critical 
value of the charge 

ge = - VH^ - 4m2) , (108) 

which separates two qualitatively different kinds of motion. In the regime of |g| < qc, 
{J^y takes both positive and negative values and r(r) obeys both tanh-type and tan- type 
solutions, representing bounded and unbounded motions respectively. Alternatively, (|108|) 
gives the critical value of H for flxed n and q, or the critical value of k for flxed H and q, 
both corresponding to the transition from bounded to unbounded motion. 

In Fig. 7 we show r(r) plots for Hq = 3 and flve different values of q. As |g| increases the 
period becomes large. When |g| exceeds the critical value q^ = 0.2700907567, the motion 
becomes unbounded and the separation of the two particles diverges at flnite r. A similar 
transition is found in Fig. 8, where r(r) plots are depicted for |g| = 0.1 and flve different 
values of Hq and the transition occurs at Hq = 7.21249. 

Before proceeding to the analysis of the phase space trajectories, it is instructive to 
compare the general structure of the determining equation for the repulsive case with that 
in no-charge case. In the q = case the determining equation (^) of the Hamiltonian for 
equal masses becomes 
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Figure 5: The phase space trajectories ol Hq = 3 for the exact, the hnear, the non-relativistic 
cases and the flat electrodynamics. 
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Figure 6: Phase space trajectories of the relativistic and the non-relativistic cases under the 
same initial condition po — "n/S/S at r = 0. 
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Figure 7: The r(T) plots for Hq — 3 and five different values of \q\ in the repulsive case. 
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Figure 8: The r(T) plots for |g| = 0.1 and five different values of Hq in the repulsive case. 
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where Y and Z are defined by 



H = Jp'^ + m^ + ep —Y{p,r), (110) 



K\r\ 



p2 _|_ ^2 _ _ (111) 



Equation ( [10 9|) has three formal solutions shown in Fig. 9: 



sol 1: line F-0 ; trivial solution Y = Z, 

sol 2: curve A-B-O-C-D; Y = W{-Z e^), Z < W-\e-^) = 0.278, 

sol 3: curve E-F-G; Y = W{Z e^), Z < 0. 

Here W{x) is the Lambert W function defined via 

y.ey = x y = W{x). (112) 

Since Z > for real p, only sol. 2 represents a physical solution. (Sol.l holds for the special 
case of massless particles with no interactions.) 

Of the four types of the determining equations (|76D -([79|) in which J7a is replaced by JT", 
tanh-type A is a generalization of sol. 2 and tanh-type B is that of sol. 3. Tan-type A and 
tan-type B are new equations appearing specifically in the repulsive case and sol.l is a special 
case of tan-type A and B solutions with JT" = and g = 0. Tan-type A and tan-type B 
trajectories yield a countably infinite series of unbounded motions of the particles. This is 
in strong contrast with the Newtonian case, in which only one trajectory exists for fixed H 
and |g|. 

Fig. 10 is a diagram of the physical region of (|(3'|,p) parameter space in the case of i^o = 3. 
The shaded area of > and S > is the region where tanh-type A and tanh-type B give 
the actual trajectories. The boundary of this area is fixed hj p = ±po = ±^ {Hq/2Y — rn? 

and p = \q\ + Hq/2. The values of |g| at the intersections of these boundary lines 

are denoted as Qc and qa, of which is the critical value ( |108| ). In this area the tanh- 
type B solution is realized in a quite narrow region between p = \q\ + Hq/2 and 
J' — B = (dashed curve). The motions of tan-type A, B are realized in the area of J'^ < 
whose boundaries are p = ±^2/k \q\ + Hq/2. Consider a |g| = const line (the dotted 
line in Fig. 10) and define by pi and p2 the momenta of the intersections of the line with 
p = \q\ + Hq/2 . 

For the case of < |g| < gc the allowed value of p is divided into two parts: —pQ < p < Pq 
and p2 < p < pi. The solution in the former region is tanh-type A and the solutions in the 
latter region are both tan-type A and B. We show in Fig. 11 the phase space trajectories 
for Hq = 3 and |g| = 0.25, in which the solid curve (A^) denotes the bounded motion given 
by tanh-type A, and the dotted {An : n = 0,1,...) and the dashed (i?„ : n = 1,2,...) 
curves represent the infinite series of unbounded motions specified by tan-type A and tan- 
type B, respectively. For the unbounded motions pi and p2 are the asymptotic values of 
the momentum and the two particles simply approach one another at some minimal value 
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Figure 10: The diagram of the physical region of for ifo = 3 in the charged repulsive 

case. 
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of \r\ and then reverse direction toward infinity. In the figure we added also the trajectory 
of flat-space electrodynamics (a dot-dashed curve), which is the only solution of the theory 
for given values of Hq and |g|. The motion of tan-type A, B has a specific feature that r(T) 
becomes infinite at a finite proper time (but an infinite coordinate time). In Appendix C we 
shall present a simple model in fiat space-time that has this feature. Fig. 12 shows the r(T) 
plots for Hq = 3 and |g| = 0.25. 

As \q\ approaches Qc the trajectories 'W " and "Aq " come close to one another, meeting 
at \q\ — Qc, and then for |g| beyond the critical value they form new unbounded trajectories 
"Nl" and 'W2" as shown in Fig. 13. This latter case will be discussed in the next subsection. 
It should be stressed that the existence of two types of motion for fixed H and g is a new 
aspect of the relativistic gravitation theory, and has no non-relativistic analogue. 

5.3 The repulsive case: the large \q\ regime {\q\ > Qc) 

For |g| larger than qc the electric repulsive force overwhelms the attractive gravitational one 
and only unbounded motion is allowed. For the case of gc < \q\ < Qd the allowed value of 
p is —po < p < pi. Fig. 14 shows the phase space trajectories of ifo = 3 and \q\ = 0.3. 
Here the new unbounded trajectories A'"l and A^2 are reahzed instead of N and ^40. The 
solution corresponding to the shaded area —po < p < p2 is N2, while the solutions for 
P2 < P < Pi are A^l, An and Bn {n = 1,2, . . .) and they are all described by tan-type A 
and B. The trajectories in r < region are obtained from those in r > by replacing the 
signs of both r and p. By comparing all these trajectories with the analogous trajectory in 
flat space electrodynamics(a dot-dashed curve), we can see how the effects of gravity deform 
the flat-space trajectory. (Remember that there exist also time-reversed trajectories with 

6=-l.) 

The phase space trajectories for the case of |g| = 2 > qa are depicted in Fig. 15. Since p2 
is smaller than —po all solutions are tan- type A and B, and a characteristic cusp appears at 
r = in the trajectories A'"l and A'"2. In the flgure the trajectory specifled by the symbol 
BOl is a combination of Bo + Bi, indicating that the solution switches between Bq and Bi, 
namely, Bq for —po < p < po and Bi for P2 < P < — Po and Po < p < Pi- Similarly S12 is 
composed of a combination of Bi + B2. 

5.4 H < 2m case for repulsive charges 

In 1 + 1 flat-space electrodynamics we know that for attractive charges the total energy of 
the system is restricted to H > 2m, but for repulsive charges no restriction on H exists. 
Unbounded motion is also realized for H < 2m (the explicit solution is presented in Appendix 
B) . It is expected that in a general relativistic theory the restriction on H is identical. ^From 
the (|q'|,p) diagram in Fig. 10 we know that for H < 2m the shaded area disappears and only 
the region of J'^ < remains for the unbounded motions. In Fig. 16 we show the phase 
space trajectories for = 1 and |g| = 1. All types of unbounded motions An and Bn are 
realized. The unbounded trajectories Nl and N2 that appeared in if > 2m turn to Aq 
again. As compared with the flat-space trajectory (a dot-dashed curve) all trajectories are 
curved more toward the r-axis (due to the additional effect of gravitational attraction) and 
are shifted in the direction of the positive p-axis. This shift is caused by the p-dependence 
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Figure 11: Phase space trajectories of the bounded and the unbounded motions for Hq — 
3, m = 1 and \q\ — 0.25. 
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Figure 12: The r(T) plots for the parameters Hq — 3,m — 1 and |g| = 0.25. 
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Figure 13: Transition from tlie bounded motion to tiie unbounded motion across \q\ = Qc — 
0.2700907567 for Hq = 3. 
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Figure 14: Phase space trajectories of the unbounded motions for Hq — 3 and |g| = 0.3. 
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Figure 15: Phase space trajectories of the unbounded motions for Hq — 3 and \q\ — 2. 
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of the gravitational potential. As remarked previously there exist corresponding trajectories 
with e = — 1 and invariance under time-reversal is retained. 



6 ELECTRODYNAMICS WITH Ae ^ 

In a recent paper P] we presented a detailed analysis of two-body motion with no charge, 
namely particle dynamics in lineal gravity. The effects of a cosmological constant (Ag = A) 
may be incorporated into a momentum-dependent potential between particles, as shown in 
the two parameters expansion in terms of k and A^/k"^ 



2 



H = 2\Jp'^ + m'^ + —{^Jp'^ + m'^ — epY \r\ + —{y p'^ + m'^ — ep)^ r"^ 

, r— 4 3 Ae ep K ep ^ 

A^ en 

+ -^- , , , ^,,3,, \r\ + ---. (113) 

The exact phase space trajectories (r(r),j9(r)) indicate that a negative cosmological con- 
stant A < acts effectively as an attractive force leading to bounded (periodic) motions, 
which are specified by the tanh-type A equation ([76|), whereas a positive cosmological con- 
stant acts effectively as a repulsive force. One noteworthy special situation takes place for a 
particular range of negative A and small m: both r(r) and the phase space trajectory have 
a double peak structure. An example is shown in Figs. 17 and 18 for Hq = 10, m = 0.02 
and A = —0.5. Two particles starting at r = with initial momenta in opposite direc- 
tions depart one another, reach a maximum separation, and then reverse direction due to 
the attractive force. However at some point they reverse direction again, reaching a second 
maximum before returning to the starting point. This peculiar behavior takes place due 
to the induced p-dependent A potential combined with the gravitational attraction and the 

2 4 

relativistic effect. For < A < Ac = 2(^2 ^m^) ' bo^^ bounded and unbounded motions are 
realized for a fixed value of H, as with the motions depicted in Fig. 11. For A > Ac only 
unbounded motion is realized analogous to the motions shown in Fig. 14. 

In the general case of electrodynamics with a cosmological constant the dynamics of 
particles is governed by a combination of four factors: gravitational attraction, the electric 
force between charges, the effect of the cosmological constant and relativistic effects. The 

solution is characterized by the signs of 7^ and (^Hq — mf{T) + j^'^ — 8^6162 — 8Ae. We 
shall investigte the motion in four combinations of the signs of Ag and the charge, separately. 



6.1 Ae < and attractive charges 

In this case all interactions (gravitational, cosmological and electric) are attractive and J'^ is 
positive. The motion becomes necessarily bounded and described by the tanh-type solution. 
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Figure 16: Phase space trajectories of the unbounded motions for Hq < 2m {m — 1, Hq ~ 1 
and \q\ — 1). 
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Figure 17: The r(r) plot for Hq = 10,m = 0.02 and A = —0.5 showing a double peak 
structure. 
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Figure 18: The phase space trajectory corresponding to the r(T) plot in Fig. 17. 
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The period of the bounded motion is 



7m > 0, 

7m = 0, (114) 



7m < 0. 



For negative Ag motion is allowed for total energy larger than \J\Ke\/ k^. As the total energy 
increases both the period and the amplitude of the motion become large and the trajectory 
r(r) deforms just as shown in Fig.l in Ag = case. 

The main purpose in this sub-section is to investigate the effects of Ag on the motion, 
especially on how the double-peak structure appears in the system of charged particles. We 
find that the double peak is caused by the interplay among the A potential, gravitational 
attraction and relativistic effects, and is suppressed as the attractive force between charges 
becomes strong. Fig.19 shows the r(T) plots for Hq = 10,m = 0.02, |g| = 0.1 and four 
different values of negative Ag. The double peak appears for Ag = —0.5. As Ag approaches 
its lower bound the form of the phase space trajectories changes from an ^'-shaped 

curve to a double peaked one and then to a diamond shape as depicted in Fig.20. In Fig.21 
we trace how the double peak in Fig.19 is affected by the value of charge We see that 
for small values of \q\ the double peak structure survives, but it disappears for large \q\. 
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tanh~"^ i 2po77;^(i+y/7F)g 

V [H{l+^)-'1e^Jpl+m??-lipl-% 



32po(l+v/7F)g 



Km 



2pov^^7;;r(l+v/7H)g 



^^tan-1 



6.2 Ag < and repulsive charges 

For the case of repulsive charges may become negative. We can classify the solutions in 
terms of the diagram as was used in the sub- section 5.2. The boundary of = is 

given by 

' ^ ' 2 / „ A~ 



± 



(115) 



2 J \\k 

There are two types of{\q\,p) diagram depending on whether the vertex {qq — ^/^, \j (f-) — 
oi Jl — Q is within the region — po <P<Po'-Po — y^(i?/2)2 — m?, or not. Fig. 2 2 is the di- 
agram for m > Y^2|Ag|/«; where the vertex is outside the region —po < p <po- The physical 
region consists of the shaded area oi Jl > and B > 0, and the area of < 0. The former 

area is the region of tanh-type solutions and the latter is the region of tan-type solutions. A 
narrow region between J'\~ B = (dashed line) and = in the shaded area is the region 
corresponding to tanh-type B solution. The Qc and are the values of |g| at the intersections 



of Jl = with p = ±po, namely, = {^(i//2)2 - 2|Ag|/«:2 _ ^(i//2)2 - m^} + |Ag|/«: 



and = {^(///2)2 - 2\A,\/k^ + ^{H/2Y - w?} + |Ag|/«:. 

The motions corresponding to this diagram are classified into three categories according to 
the \q\ value: 
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Figure 19: The r(T) plots for Hq — 10,m — 0.02, |g| = 0.1 and four different values of Ag. 
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Figure 20: Phase space trajectories corresponding to the plots in Fig. 19. 
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Figure 22: The diagram for Ag < and repulsive charges under the condition m > 

^J2\h.e\/ K. The dashed hne represents J^ — B — 
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(i) IqI ^ Qo '■ bounded motion, 

(ii) go < \q\ < Qc '■ both bounded and unbounded motions, 

(iii) Qc < \q\ : unbounded motion. 



For a small |g| in the category (i), the attractive effect of Ag < is stronger than the repulsive 
effect between charges (|Ae| > nq^) and Jl is positive. The phase space trajectories resemble 
the trajectories in Fig. 4, but they become more S'-shaped as |g| increases. For the categories 
(ii) and (iii) the (Iq"!,!?) diagram is nearly the same as the diagram in Fig.lO. The phase 
space trajectories for (ii) are just like those in Fig.ll. The trajectories of the motion for 
the category (iii) are divided into two cases of Qc < \q\ < Qd and Q'd < In the case of 
?c < \q\ ^ Qd the trajectories are analogous to those in Fig. 14 in which the N2 trajectory 
is a tanh-type, while in the case of qa < \q\ all trajectories are tan-type and like those in 
Fig. 15 where A^^i and N2 have cusps at r = 0. In Fig.23 the r(r) plot for each category 
is depicted for the parameters Hq = 3,m = 1.2 and Ag = —0.1 with qo = 0.316 and 
qc = 0.491: \q\ = 0.1,0.4 and 0.6 for the categories (i), (ii) and (iii), respectively. As \q\ 
becomes large, both the period and the amplitude of the motion increase and finally the 
motion becomes unbounded, because the repulsive force between charges prevails over the 
attractive forces of gravity and the cosmological constant. Fig.24 shows the corresponding 
phase space trajectories. 

Since the analysis in sub-section 6.1 indicates that the double peak structure appears 
for negative Ag, small m, and small attractive \q\, it can be inferred that for the repulsive 
charges, if |g| is sufficiently small, the double peak still survives. The smallness of m corre- 



sponds to the condition m < y2|Ae|//t for which the (|?|,p) diagram is given in Fig. 25 and 
the vertex of = is within the region —po <p <po- The motions are classified into two 

categories: 

(i) \q\ ^ Qo '■ bounded motion or unbounded motion, 

(ii) qo < \q\ unbounded motion. 

The physical region of the category (i) belongs to the shaded area of Jl > and B > 0. 

We can find the double peak structure for a certain range of the parameters as shown 
in Fig.26, in which Hq = 10,m = 0.02, A^ = -0.5 and \q\ = 0, 0.1, 0.14. The above 



parameters of m, Ag and |g| correpond to 7^ < 0. As |g| exceeds qm = y i^^y'^l^el — K-rn'^/'2, 
the solution to jm = 0, we encounter a new situation. The maximum turning point of r 
extends to infinity and the trajectory splits into two nonperiodic motions as shown in Fig. 27 
and Fig. 28. 

For \q\ — qm the asymptotic value of p is m^{H'^ — |Ag|)/8|Ag| and in the (|?|,p) diagram 
of Fig.25 this corresponds to the vertex point of a dashed-curve oi Jp^ — B = 0. For qm < \q\ < 
qo the two asymptotic values of p are the intersections of |g| = const line with J7a — B = 
curve. The region between these asymptotic values originally belongs to tanh-type B solution 
but no trajectory exists because J^a> H — \mf — m/ f\ for the parameters in this region. 

The trajectories of the motion for category (ii) are diveded into three cases of go < |?| < 
Qc, Qc ^ \q\ < Qd and qa < \q\. For go < \q\ < Qc under the condition m < ^Jq\K]\/k , 
a, < Q region is sandwiched by the shaded areas, in contrast to m > ^Jq\KJ\/k case 
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Figure 24: The phase space trajectories corresponding to the plots in Fig. 2 3 



49 



p 




Figure 25: The diagram for Ag < and repulsive charges under the condition m < 

^J2\h.e\/ K. The dashed hne represents J^ — B — 
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Figure 26: The r > part of the phase space trajectories for Hq — 10, m = 0.02, Ag = —0.5 
and three different values of repulsive \q\. 
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|q|=qm (0.140712) 
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Figure 27: The r > part of the phase space trajectories for Hq — I0,m — 0.02, Ag 
and Qm < \q\ < Qq. 
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where the region for < \q\ < qc is separated into two parts, namely, a shaded area and a 
Jll < area. A typical phase space trajectory for go < \q\ < Qc is shown in Fig. 29 where the 
parameters are Hq — 3,m — 0.6, Ag = —1 and \q\ = 1.1. Here between two split nonperiodic 
trajectories (solid curves) there appear an infinite series of tan-type A, B solutions of the 
unbounded motions. For the motions of the cases qc < \q\ < qa and qd < \q\ the outlines of 
the unbounded trajectories are easily inferred from Figs. 14, 15 and 29. 

6.3 Ae > and attractive charges 

For the attractive charges the boundary curve ol Jl — Q is 




(116) 



which opens out to the direction of p-axis. Fig. 30 is the (Iq'l,^) diagram for 

H < ^ K^m? /2Kf, + Am?. The motions for this diagram are classified into two categories: 

(i) < |g| < go : both bounded and unbounded motions, 

(ii) go < kl • bounded motion. 

The phase space trajectories for category (i) are just like the trajectories in Fig.ll, while for 
category (ii) the motions become simply bounded as the attractive charges overwhelm the 
repulsive efi'ect of the efi^ective cosmo logical constant. 

Under the condition of H > ^ K^m? /2Ke + 4m^ the situation is slightly different from 
the previous case, as shown in the (|g|,p) diagram of Fig. 31. In this case the motions are 
classified into three categories according to the |g| value: 

(i) < |g| < gc : unbounded motion, 

(ii) gc < |g| < go : both bounded and unbounded motions, 

(iii) go < |g| : bounded motion, 

where g^ = ^Ae/re - (k/2) {^{H/2f + 2k J k"^ - ^{H/2f - rJf . Since this (|g|,p) dia- 
gram is like the inverse of the diagram in Fig. 22, categories (i), (ii) and (iii) correspond to 
the previous categories (iii), (ii) and (i) of m > ^J2\KJ\/ n in sub-section 6.2, respectively. 
In Fig. 32 the rij) plot for each category is drawn for the parameters i^o = 3,m = 1 and 
Ag = 1 with go = 1 and Qc — 0.745: |g| = 0.5,0.8 and 1.5 for the category (i), (ii) and 
(iii), respectively. As |g| increases, the motion changes from unbounded to bounded and 
then both the period and the amplitude decrease.The phase space trajectory for each r(T) 
is easily inferred from Fig. 24. 

6.4 Ae > and repulsive charges 

In this case, since the cosmological repulsion is commensurate with electromagnetic repulsion, 
classification of the motion is analogous to that of repulsive charges with A,. = in sub- 
section 5.2. Under the condition H < H?im? /2h.(. + An? the (|g|,p) diagram is given by 
Fig.33 which has the same pattern as Fig. 10. The physical region is classified into two 
categories: 
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Figure 29: The phase space trajectories for Hq — 10, m — 0.02, Ag = —0.5 and go < \q\ < Qc- 
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Figure 30: The diagram for Ag > and attractive charges under the condition 

H < ^K2^V2Ae + 4m2. 
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Figure 31: The diagram for Ag > and attractive charges under the condition 

H > K^m? /2h.e + Am?. The dashed hne represents Ja — B — 
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(i) |g| < qc : both bounded and unbounded motions, 

(ii) qc< \q\ '■ unbounded motion. 

The trajectories of the motion for the category (ii) are divided into two cases of Qc < \q\ < Qd 
and qd < \q\. The phase space trajectories and r(r) plots for |g| < qc are analogous to those 

in Fig. 11 and 12. For a larger H of H > K^m? /2Ke + 'im? the motion is always unbounded 
as inferred from the diagram of Fig. 34. 



7 MOTION OF UNEQUAL MASSES 

What time variable is adequate for the analysis of the motion of the two particles with 
unequal masses? The proper time (|^) of each particle is different as 



dTi 



dT2 

where K = = and Y = Yj^ 
case [§], we propose using 



dt 



dt 



lQYKoK2m2 



(117) 



JKM2^Jp'^ + ml 
yi. As with our analysis of the electrically neutral 



dr = dt 



IQYK< 







KiK2mim2 



1/2 



JK \MiM2Jp^ + mf 7p2 + ml 



:ii8) 



which is symmetric with respect to 1 <-> 2 and reduces to the proper time ( ^If ) when mi = m2. 
In terms of this variable the canonical equations are expressed as 



dp 
dzi 



dr 
df 



1 / K1K2M1M2 + rnl + rnl 
Ak \ mim2 



1/2 



(119) 



MiM2\Jp^ + mi^p^ + mi 
KiK2mim2 



1/2 



tJ IQ 




Y 
K 



(120) 



MiM2^p^ + ml^p 
KiK2mim2 



1/2 




Y 
K 



+ 



K2 




(121) 



Note that r still describes the proper distance between the particles at any fixed instant. 
Unlike the equal mass case, the integration / dp{KiK2MiM2\J p^ + mlJ p^ + )^^/'^ 



can not be performed within the framework of elementary calculus. We resort to numerical 
calculation for solving above equations. 



59 



p 




Figure 33: The diagram for Ag > and repulsive charges under the condition H < 

K^m? /2h.e + im?. The dashed hne represents Ja — B — 
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Figure 34: The diagram for Ag > and repulsive charges under the condition H > 

K^m? /2h.e + im?. The dashed hne represents Ja — B — 
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As in Sec. VI we analyze the motions by plotting r(r) in four combinations of the signs of 
Ae and the charges. In the case of a negative cosmological constant and attractive charges, 
the r(r) plots are shown in Fig.35 for various mass ratios mi/m2 in the fixed = 10, 1712 = 
l,Ae = —1 and |g| = 1. When the mass ratio gets larger than unity (the case of mi = 
7712 = 1); the gravitational attraction is stronger and the proper distance between particles 
as well as the period becomes shorter. As the mass ratio gets smaller gravity becomes weak; 
however for quite small mass ratios the attractive effect due to the cosmological constant 
prevails and the period changes to being short again. The effect of the cosmological constant 
is seen by changing |g| small but preserving the values of other parameters. Fig. 36 shows 
the r(r) plots for |g| = 0.1. While both the proper distance and the period becomes large 
as a whole, the double peak structures appear for a small mass ratio. 

In the case of Ag < and repulsive charges unbounded trajectories appear as |g| increases. 
An example is shown in Fig.37, in which r(r) for a small mass ratio mi/m2 = 0.1 becomes 
unbounded because the repulsive force between charges exceeds the weak gravity and the A 
force. 

In the case of Ag > and attractive charges, as we analyzed in Sec. VI the trajectory 
changes from unbounded to bounded as |g| increases. In Fig. 38 we plot r(r) for various 
mi/m2 in the fixed Hq = 10, m2 = 1, Ag = 0.1 and |g| = 0.1, from which the strong repulsive 
effect of a positive Ag is perceived. The motions for the case of Ag > and the repulsive 
charges are mostly unbounded except for the case of the large masses. 



8 STATIC BALANCE IN (1+1) DIMENSIONS 



In this section we treat the problem of static balance in (1+1) dimensions. This problem 
originated in attempts to find the exact solution of the N{N > 2) body system in Einstein's 
theory of general relativity |]12|. In (3 + 1) dimensions gravitational radiation carries away 
energy from a system of particles, and so exact solutions have so far been unobtainable. One 
simple way to search an exact solution for > 2 is to balance the gravitational attraction 
with some extra repulsive force, a natural candidate of which is the electric force. The first 
trial was achieved by Majumdar ||TB[ and Papapetrou [T^\ for N = 2 and afterwards it was 
generalized to N bodies on a line [|1^. Their condition for balance is 



rrii 



:i22) 



and is much more strict than the corresponding condition in Newtonian theory 



Gmi 777.2 ~ 



6162 



(123) 



Since then, the question has long been raised as to why the balance condition differs in 
the relativistic and non-relativistic cases. Some people conjectured that there should 
be an exact solution in general relativity under the condition ( |123| ). Others |T^, |TB| showed 
in the (2nd) post-Newtonian approximation that the condition ( |123| ) is incompatible with 



the static balance condition in general relativity. Clearly (|122D is a sufficient condition for 
static balance but there exists no proof that ( |122|) is a necessary condition. On the other 
hand a test particle analysis [0 suggested that the condition (|123|) was neither necessary 
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Figure 37: r(T) plots for different values of the mass ratio raxjm^ for = 10, m2 = 1, Ag = 
— 1 and ei — 62 — \q\ — 0.6. 
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Figure 38: r(T) plots for the different values of the mass ratio mi/m2 for Hq = 10, m2 = 
1, Ae = 0.1 and ei = —62 = \q\ — 0.1. 
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nor sufficient, but a separation-dependent balance position might exist in general relativity. 
Recently several numerical trials on the separation-dependence were reported but 
so far no one has been able -in any relativistic theory of gravity- to succeed in finding 
analytically an exact solution under ( |123| ) or another (separation-dependent) condition. 

In (1+1) dimensions the absence of radiation makes it easy to fix the balance condition 
in terms of the determining equation (0). From the relation 

or kJ 
the balance condition dH/ dr = leads to 



K1K2 = {2Ko + 2K- i^^Jp^ + mj) {2Ko + 2K ~ i^^Jp^ + m^) = , (125) 

where K = K± = kX and H = -K. In deriving (|125|) we excluded the possibility that 
Kq = because it gives an unphysical solution. 

It is evident from (^) that the condition (|125| ) means at the same time M1M2 = 0, 
namely, 



{2Ko -2K + K,Jp^ + mj) {2Ko -2K + n^p^ + mi) = , (126) 
The equations ( |125| ) and ( |126| ) lead to 



+ {2K - K\jp^ + m\){2K - K\jp^ + m^) = (127) 
and 



Kq{^K - k\Ip^ + mf - k\Ip^ + mi) = . (128) 

The solution to (|12|) is 

E = ^pi + ml + y^p2 + ^2 (^29) 
and the insertion of this solution into ( |127| ) leads to 



2 ^ y p2 _,_ ^2 _ ep){^Jp'^ + mi - ep) - 6163 = . (130) 

This condition is the force-balance condition and fixes the value of momentum p = Pc = 
const. Under the condition (|130|) the two particles move with a constant velocity. The 
condition ( |130| ) and the Hamiltonian ( |129| ) indicate that the full Hamiltonian must have the 
simple structure 

H = \Jp'^ + mf + y^p^ _^ Tpfi'^j^i^-i^^pi -\- mj- ep){\Jp^ + mi - ep) - 6162 1 F(|r|,p) . (131) 

Actually, in the case of Ag = the K-expansion of the Hamiltonian leads to 

1 



H = y p'^ + mj + y p'^ + ml - -6162 \r\ 



+^ I V^P^ + "^1 + "^i \r\ — ^pi^/p"^ + ml + ^jp^ + mi )|r| + p^|r| 



-eie2{\lp'' + m{ + yp^ + m^2)r + -eeiCsP r +^(6162) \r\ 
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K ( 1 I / 1 

^42 12^^^*'^^^ + + ^p'^ + ml- 2ep - -eie2\r\f 
+—^—6162 \r\^{\/p'^ + mj + Jp'^ + ml - 2ep - l-eie2\r\f 



1 2 

r 



6 X 42 



12( Jp2 + rnl- Jp^ + miy + (eiea)^^ 



X (yp2 + 77^2 + y'p2 _^ ^2 _ 2ep - -eie2|r|) 



+ 



3 X 42 

From ( |132| ) we get explicitly 



^eie2 \r\%^p^ + ml - ^p^ + m^ f - -J-—{e,e2)^\r\'] + 0{k') (132) 
< 4^ V V 15 X 4^^ J 



F(|r|,p) = -|r| + -(y'p2 + rnl + ^p^ ^m\- 2ep) + 0{\r\^) . (133) 
The solution to the condition ( |130| ) exists only for 6162 > 0; 



(t) 



2 

22 221 



Pc = ± ■ (134) 

v/^Ke^J (f mf + eie2)(f + 6162) 



When the particles are initially at rest {pc = 0), the condition ( [1301 ) becomes 

^mim2 - 6162 = . (135) 

This is the condition of static balance and is identical with the condition in Newtonian theory 
in (1+1) dimensions. Note that in Newtonian theory ( |135| ) is the force-balance condition 
which includes both the static case and a uniform motion. However in the relativistic case 
(|135|) represents only the static balance condition - the condition of force-balance (|130|) in 
general depends on the momentum and allows a uniform motion in the C.I. ststem in which 
two particles approach or recede with the designated constant momentum ( |134| ). 

The above result may seem to provide suggestive evidence that in (3+1) dimensions (|123|) 
could be a sufficient solution. However, the post- Newtonian Hamiltonian for the system of 
two charged bodies in (3+1) dimensions is given as 

Pi , pI {pI? {pI? 



H = mi + 1712 H 1- q 

2mi 2m2 8mf 8m2 

Gmim2 ^^^3/p^^pl\ 7 (p^ ■ Pg) {Pi-r){P2-r) 



r [ 2 ymf ml J 2 mim2 2mim2r'^ 
^ 6162 L_ {Pi ■ P2) _ jpi ■ r){p2 ■ r) \ ^ G'2mim2(mi + 7722) 



Anr y 2mim2 2mim2r'^ J 2r2 

G(mi+ 7712)6162 , G(?Tii6| + 77126?) 

~A — 2 + Q — 2 ' ^^^^> 

which is derived from Bazahski Lagrangian and also the ADM canonical formalism. The 
structure of the Hamiltonian (|136|) is rather complicated compared to (|131|) . For example 
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in (1+1) dimensional theory the charges appear always in a combination 6162, while in 
(3+1) dimensions there is a combination of mie| + m2e\ and in higher approximations more 
complicated combinations of mass and charge appear It is clear that the Hamiltonian 
( |136| ) supports the balance condition ( 122 ) for the static case, while it does not correspond 
to a uniform motion. What causes this difference between (1+1) and (3+1) dimensional 
theories? Is it intrinsic to the dimensionality, or is there still a possibility for satisfying 
the balance condition ( [123| ) in (3+1) dimensions? Or is there a uniform motion in the C.I. 
system in (3+1) dimensions? These are interesting open problems for further investigation. 



9 CONCLUSIONS 

Since in (1+1) dimensions the degrees of freedom of both gravitational and electromagnetic 
radiation are frozen, one expects the motion of a set of charged particles in curved- 
spacetime with a cosmological constant to be described by a conservative Hamiltonian. And 
this is what we find to be the case. We began by canonically reducing the charged A^- 
body action to first-order form and then for a system of two charged particles derived the 
determining equation of the Hamiltonian from the matching conditions. The canonical 
equations of motion given by the Hamiltonian can be solved exactly when they are expressed 
in terms of the proper time, and we have given the explicit solution. To our knowledge this 
is the first non-perturbative relativistic solution of this problem. 
We recapitulate the main results of this paper: 

(1) In (1+1) dimensions the square of the electric field plays the same role as the consmo- 
logical constant and an overall constant part is incorporated into the effective cosmological 
constant Ae = Aq — /^(Z^a ea)^/4, which induces the momentum-dependent potential in the 
Hamiltonian. Effectively Ae > acts on the particles as a repulsive potential and Ag < 
acts as an attractive potential. 

(2) The theory in the case Ag = is a general relativistic electrodynamics which is an exten- 
sion of fiat space Newtonian electrodynamics. For the attractive charges the motion becomes 
bounded similar to |g| = case with smaller amplitude and period. For repulsive charges 
the electric force between two particles competes with the gravitational attractive force. For 
|g| < Qc and a fixed energy not only bounded motion but also an infinite series of unbounded 
motions is realized. Such multiple solutions do not exist in Newtonian theory. For |g| > 
all motions become unbounded. 

(3) The condition for static balance in (1+1) dimensional general relativistic electrodynam- 
ics is identical with the Newtonian balance condition in the fiat space. This simple result 
contrasts strikingly with the unresolved situation in (3+1) dimensions. 

(4) In spacetimes with nonzero Ag the motions are more complicated. For a fixed total 
energy Hq and Ag the solutions are classified in terms of (|(3'|,p) diagram in which the region 
of the parameters for the bounded and/or unbounded motions are depicted in terms of the 
curves of p = ipo, J]i = ^ and J7a — -B = 0. By drawing a g = const line on the diagram we 
can easily find what types of the motion are realized. 

(5) For a certain range of Ag < and a small |g| the double peak structure appears in r(r) and 
the phase space trajectory, which is caused by a subtle interplay amongst the momentum- 
dependent A potential, the gravitational and the electric potentials, and relativistic effects. 
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(6) For unbounded motion the mutual separation of two particles diverges at finite proper 
time. This is common for the cases of Ag = and Ag 7^ 0. 

(7) In the unequal mass case the basic features of the motion are the same, but a double 
peak structure appears more clearly than in the equal mass case. 

A number of interesting questions arise from this work. First, it would be of interest 
to know how the features of the motion we have found are modified as one increases the 
number of bodies in the system. For neutral bodies this situation is relevant in modelling 
stellar systems and galactic evolution, where non-relativsitic one-dimensional self-gravitating 
systems are employed as tool toward understanding such dynamics. For charged bodies the 
relevance is most likely in terms of the physics of the early universe, where charged black 
holes, strings and domain walls interact in a highly relativistic setting. Second, it would 
be interesting to study the circumstances under which black hole event horizons can form - 
we expect this will involve investigating the regime where eq. ( PP| ) is not satisfied. A third 
possibility would involve generalizing our work to other (1+1) dilatonic theories of gravity. 
Moreover, a full quantum treatment of the problem would also be of considerable interest. 
We intend to turn our attention to these problems in the future. 



APPENDIX A: SOLUTION OF THE METRIC TEN- 
SOR 

Under the coordinate conditions (|35D the field equations (pT]), (pH]), and (|32|) become 



3k 



7r + Noi—7r' + —(^0 + - V - ^ 



A, 



K 



2Jpl + ml 



2 J pi + mi 



--5{x - Z2{t)) 



+Ni {tt' + pi5{x - zi{t)) + p25{x - Z2{t))} + — A^o^' + A^^tt = 
kttNo + A^{ = , 

^ + 2kNott - Ni^' = . 



(137) 
(138) 
(139) 
(140) 



The procedure to solve these equations is just the same as that in the case of \^ = 0, which 
is given in the Appendix of the previous paper 0. 

The solutions have the same form as those of V = when they are expressed in terms of 
K±, Kq, Ki^2,-Mi,2, y± and J. For r = zi — Z2 > the metric tensor is determined as 
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8(Y^ Y_\KoK, 
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2KoJ 
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8 /r+ , r_\ KoKi 



K+ 1 J \K+ ' K^l Ml 



+ — I e 

J \A + 

/ Yo (Y^ Y. 



A"+(a;-2i) 



- 1 
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Y^ 



Y^ 



1 



Ko K_ / 



K1K2 - -^KiK2{zi + Z2) \ . (142) 



The metric tensor for r < is simply obtained by interchanging the suffix 1 and 2. With 
these metric tensors and the canonical equations, the field equations ( |137| ), ( |138| ) and ( |139| ) 
are proved to hold in a whole x space. 

As we showed in the previous paper, to satisfy ( |14CI| ) the dilaton field \l/ needs an extra 
x-independent function /(t), which has no effect on the dynamics of particles. After lengthy 
calculation Eq. (|14CI|) leads to 
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--iKoiz,-Ko2Z2) + -[- 
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2KoKi 
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Ml M2 



(143) 



Thus fit) is uniquely determined. 



APPENDIX B: LINEAR APPROXIMATION OF THE 
EXACT SOLUTION IN Ae = 

In this appendix we investigate how the exact solution is related to its corresponding solution 
in Newtonian theory. Take the tanh-type A solution, namely a periodic motion with 6162 < 
or small 6162 > 0. For motion starting from r = at r = tq = 0, the relationship between 



the energy and the initial momentum is Hq = 2^Jpq + m?. 
In the n expansion, feij) is expressed as 



fe{r) 



Pq + m? + epo 



e2m 



£162 



—Km 



Pq + 
6162 



(e^'''^ -l)-{^pl + ^^ + epo\ 



er 
Am 



6 2™ 



+ Oin^llU) 
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and p{t) is 



P 



eie2T 



Pq + sinh h po cosh 



2m 



eie2T 
2m 



Km 



2e 



e ™ 



6162 



X 



pI + m^ 



6162 



1 - 



pl + m^ + epo) T— e 2" 
4m 



+ .(145) 



The trajectory r(r) > in the hnear approximation becomes 



r r 



6162 
6162 



+ ( 1 — cosh 



2m 



Po sinh 



eie2r 
2m , 



Po + smh h Po cosh 

2m 2m 



Po + rn?{Jpl + m2 - epo) 



4fi; 



6162 

X 



Po + m^ ( 1 — cosh 



6162 
2m 



— Po sinh 



1—62' 



-6162 



emr 



2m 



2eie2 



Po + m"* 



e sinh 



eie2T 
2m 



epo cosh 



eie2r 



2m 



1 



66162 

Furthermore in the hmit k — > the solution is 



^f. v,^i^2r\ . , ei62rl2~ 

Po + m^ 1 — cosh— — I — Po smh— — 



2m / 



2m 



(146) 



p 



Po + m^ smh h po cosh 

2m 2m 



(147) 



6162 



Po + m2 1 — cosh 



2m 



Po sinh 



6i62r 



2m 



We can transform it into an expression in terms of the original time coordinate t, by using 
the relation 



t 



2^p^ + m2 _ ei62r 2po / , 6i62r 
— smh 1 cosh 1 



6162 



2m 



6162 V 



2m 



which is obtained by integrating (^J). Then 

6162 



p(t) 
r{t) 



Po + 



-t 



6162 



|y(Po + + m2 - y^p2 + m2| 



(148) 

(149) 
(150) 



This solution is identical with that derived from the Hamiltonian for flat-space electrody- 
namics in (1+1) dimensions: H = 2\/p^ + m2 — eie2|r|/2. The same solution can be applied 
to the system of arbitrary charges as far as Hq > 2m. 
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On the other hand, for the repulsive charges (6162 > 0) there exists the solution for 
Hq < 2m, which denotes an unbounded motion and is given by 



Pit) 
r{t) 



6162 

i 

2 

4 

6162 



(6162)2 Ho 
— — + m^- — 



(151) 
(152) 



This solution is also derived from the tan-type A solution in the limit n ^ 0. 



Another limit of the solution ( |146|) is to take Cq — 0, which leads to 



r[T) = T — K 

m 



Pq + TIT? 



T 



(153) 



In the non-relativistic approximation (keeping the terms to the lowest order ofpo/m) we get 
the solution for the motion in Newtonian gravity in (1+1) dimensions: 



r(t) = — t 

^ ' m 4 



(154) 



APPENDIX C: CAUSAL RELATIONSHIPS BETWEEN 
PARTICLES 

One of the striking features we found in the analysis of the two body motion is that in 
the repulsive trajectories Ni,N2, An and Bn the particles reach the asymptotic regime (r — 
±00, p = finite) at some finite proper time Toq. For example for A^^i trajectory it is 



To, 



Km^'Jrn 



:log 



H{i + ^/t;^) + |p+ + ^pl + m^} (^/t;;; - 7e) 



(155) 



with p+ = ^ {Vti'^H-' + SAe + V8Ae + 86162}. 

We can understand this feature in a simple flat-space model. Consider the 2-velocity 



i>^={f{aT),^Piar)-l) with dr' 
where ficrr) is some function. This 2-velocity means 



dt — dx 



dx 
Itt 



f 



and leads to the 2-acceleration 



duf" 
dr 



/V(l, 



/ 



(156) 



(157) 



:i58i 



where /' = df{T)/dr. There exist functions / such that / — > 00 as r ^ Tqc, t^o being finite; 
the particle thus becomes lightlike in a finite amount of proper time, but an infinite amount 
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of coordinate time t = J^^ drflar). An example would be / = sec (err). The acceleration is 
not constant, but increases as a function of proper time, diverging at r = Too = This is 
the situation we encounter for the unbounded motions of two charged particles (or of two 
neutral particles with sufficiently large positive cosmological constant). 

To see the causal relationships between particles and the space-time structure we try to 
pursue the path of light signal emitted from particle 2 in the metric given in Appendix A. 
The path x{t) of the light is governed by dr = 0, which reads 

The equation of the hght signal emitted inward (directed to particle 1) is 

dx 

— = No{x{t),zi{t),Z2{t),p{t))-Ni{x{t),zi{t),Z2{t),p{t)) , (160) 

and the light emitted outward is described by 
dx 

— = -No{x{t),zi{t),Z2{t),p{t)) - m{xit),z,{t),Z2{t),p{t)) . (161) 

Let's take first a typical bounded motion in the case of Ag = and the attractive charges 
for the parameters of ifo = 3,m = 1 and q = 1- In Fig.39 the trajectories of light signals 
emitted from particle 2 at various times T are plotted. The two particles are always in 
causal contact, because the inward light signal from particle 2 reaches particle 1. We see one 
striking feature for the case of Ag = 0, namely that the outward light pulse has a constant 
velocity c, since the equations (|160D and (|161|) become dx/dt = ±1, where the plus and the 
minus signs correspond to the light emitted outward from the particle 1 and 2, respectively. 

Next we look into the effect of the cosmological constant on the path of the light signal. 
Fig. 40 shows the trajecories of light for the case of a negative cosmological constant Ag = — 1 
and the same values of the parameters HQ,m and q. The outward light behaves as if it is 
subjected to a repulsive force, in significant contrast to the particles' trajectories on which 
the cosmological constant acts effectively as an attractive force as we analyzed in Sec. VI. On 
the other hand, while the positive cosmological constant causes a repulsive force between 
particles, the outward light behaves as if it undergoes an attractive force and approaches a 
kind of horizon given by the line No{x, zi{t), Z2(t),p{t)) ± ^"1(0;, zi{t), Z2it),p{t)) = which 
is shown in Fig. 41 as a narrow solid line for the case of Ag = 0.5. 

When the cosmological constant exceeds a critical value, the particles' motion becomes 
unbounded and the light signals emitted from the particle 2 exhibit new characteristics as 
shown in Fig. 42 for the case of Ag = 2.5. For small T(T < 0.9), the particles are in causal 
contact (the inward dotted curve), but for T ^ 0.9 the signal just barely catches up with 
particle 1, which is almost in light-like motion (the inward dashed curve). For T = 2 the 
inward world line x{t) is parallel to Zi{t) at large t and in the outward direction it goes 
nearly on the same trajectory with the particle 2. For large T (T > 2) the particles are out 
of causal contact with each other : a light ray sent from particle 2 toward particle 1 receives 
a strong repulsive effect and ultimately reverses direction, following behind particle 2. 
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Figure 39: The trajectories of light signals emitted from the particle 2 for the case of Ag = 0. 
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Figure 40: The trajectories of light signals emitted from the particle 2 for the case of Ae = — 1. 
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Figure 41: The trajectories of light signals emitted from the particle 2 for the case of Ag = 0.5. 
The narrow solid line denotes No{x, Zi{t), Z2{t),p{t)) + Ni{x, Zi{t), Z2{t),p{t)) — 0. 
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Figure 42: The trajectories of light signals emitted from the particle 2 for the case of Ae = 2.5. 
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